Multiple pairwise comparisons
# Load results of multiple pairwise comparisons from 12_Organ_Bias.Rmd
## Mean expression
comparisons.mean.expr.list <- vector(mode="list", length=3)
names(comparisons.mean.expr.list) <- c("observed", paste0("rnd",params$seed.run1), paste0("rnd",params$seed.run2))
comparisons.mean.expr.list[["observed"]] <- readRDS(
file.path(params$input.rds.path,
paste(params$species, "comparisons.mean.expr.rds", sep="_")))
comparisons.mean.expr.list[[paste0("rnd",params$seed.run1)]] <- readRDS(
file.path(params$input.rds.path,
paste(params$species, "comparisons.mean.expr", paste0("rnd", params$seed.run1,".rds"), sep="_")))
comparisons.mean.expr.list[[paste0("rnd",params$seed.run2)]] <- readRDS(
file.path(params$input.rds.path,
paste(params$species, "comparisons.mean.expr", paste0("rnd", params$seed.run2,".rds"), sep="_")))
# Log-coefficient of variation
comparisons.log2cv.list <- vector(mode="list", length=3)
names(comparisons.log2cv.list) <- c("observed", paste0("rnd",params$seed.run1), paste0("rnd",params$seed.run2))
comparisons.log2cv.list[["observed"]] <- readRDS(
file.path(params$input.rds.path,
paste(params$species, "comparisons.log2cv.rds", sep="_")))
comparisons.log2cv.list[[paste0("rnd",params$seed.run1)]] <- readRDS(
file.path(params$input.rds.path,
paste(params$species, "comparisons.log2cv", paste0("rnd", params$seed.run1,".rds"), sep="_")))
comparisons.log2cv.list[[paste0("rnd",params$seed.run2)]] <- readRDS(
file.path(params$input.rds.path,
paste(params$species, "comparisons.log2cv", paste0("rnd", params$seed.run2,".rds"), sep="_")))
# Residual variation
comparisons.resid.var.list <- vector(mode="list", length=3)
names(comparisons.resid.var.list) <- c("observed", paste0("rnd",params$seed.run1), paste0("rnd",params$seed.run2))
comparisons.resid.var.list[["observed"]] <- readRDS(
file.path(params$input.rds.path,
paste(params$species, "comparisons.resid.var.rds", sep="_")))
comparisons.resid.var.list[[paste0("rnd",params$seed.run1)]] <- readRDS(
file.path(params$input.rds.path,
paste(params$species, "comparisons.resid.var", paste0("rnd", params$seed.run1,".rds"), sep="_")))
comparisons.resid.var.list[[paste0("rnd",params$seed.run2)]] <- readRDS(
file.path(params$input.rds.path,
paste(params$species, "comparisons.resid.var", paste0("rnd", params$seed.run2,".rds"), sep="_")))
# Variability ranks
comparisons.expr.var.list <- vector(mode="list", length=3)
names(comparisons.expr.var.list) <- c("observed", paste0("rnd",params$seed.run1), paste0("rnd",params$seed.run2))
comparisons.expr.var.list[["observed"]] <- readRDS(
file.path(params$input.rds.path,
paste(params$species, "comparisons.expr.var.rds", sep="_")))
comparisons.expr.var.list[[paste0("rnd",params$seed.run1)]] <- readRDS(
file.path(params$input.rds.path,
paste(params$species, "comparisons.expr.var", paste0("rnd", params$seed.run1,".rds"), sep="_")))
comparisons.expr.var.list[[paste0("rnd",params$seed.run2)]] <- readRDS(
file.path(params$input.rds.path,
paste(params$species, "comparisons.expr.var", paste0("rnd", params$seed.run2,".rds"), sep="_")))